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Abstract. We calculate the S-wave scattering lengths for charmed mesons scattering off Goldstone bosons and explore 
their quark mass dependence using chiral perturbation theory up to next-to-leading order as well as a unitarized version 
of it. The quark mass dependence of all scattering lengths determined in a recent lattice calculation can be reproduced 
by the unitarized version. We also discuss signals of possible bound states in these observables. 
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1 Introduction 

The interaction strengths amongst hadrons are fundamental quan- 
tities of the strong interactions. There are two different ap- 
proaches to analyze the low-energy scattering amplitudes in 
a model-independent manner. First, one can use the chiral per- 
turbation theory (ChPT) representation of the pertinent scatter- 
ing amplitude which is given in terms of a few LECs. These 
might either be determined directly from experiment or from 
full QCD calculations performed on a discretized space-time 
lattice. A beautiful example in this context is the accurate de- 
termination of the pionic LECs I3 and from lattice simula- 
tions — see Refs. [1] and references therein. The second ap- 
proach to address the issue of scattering lengths is the exploita- 
tion of Liischer's formula I12I3I in the analysis of the lattice 
gauge theory results. It allows one to relate two-particle scat- 
tering lengths to the behaviour of certain energy levels in a fi- 
nite volume. In full lattice QCD, the maximally stretched states 
for TTTT, ttK and KK scattering have been systematically and 
accurately analyzed by the NPLQCD collaboration based on an 
extension of the Liischer formula (for a review, see Ref. |j4j and 
references therein). 

In this paper we will discuss another sector, where the same 
program can be carried out, namely the scattering of Gold- 
stone bosons off _D-mesons. A recent lattice study for these 
systems is given in Ref. |5|. However, contrary to the tttt sys- 
tem, here an additional complication may occur, for the interac- 
tion in some channels might be sufficiently strong to show non- 
perturbative phenomena even on the level of hadron-hadron 
scattering, eventually leading to hadron bound states — the 



so-called hadronic molecules. A famous example of relevance 
here is the D*g (2317) identified as a DK bound state in various 
studies 06I7I8I9I . We will discuss to what extent non-perturbative 
phenomena might already be present in the channels inves- 
tigated on the lattice and also point out the most interesting 
channels not yet studied. Note also that a detailed knowledge 
of Goldstone boson-D-meson scattering is not only of inter- 
est on its own, it is also an essential input for calculations of 
D-mesons in matter 1 10|. 

Preceding the direct lattice calculations of the Goldstone 
boson-D-meson scattering lengths, the one for the 1 — 1/2 
Dtt channel (where / denotes the total isospin) was extracted 
using the lattice calculations of the Dtt scalar form factor in 
semileptonic D n decays in Ref. [11]. The result was A/^a = 
0.29(4), i.e. a = 0.41(6) fm (where a denotes the scattering 
length). Using the same method, the value of the isoscalar DK 
scattering length was found to be effectively infinite, which was 
interpreted as a signal of a bound state of DK at threshold. 

In this paper, we present a systematic study of the S-wave 
scattering lengths of the Goldstone boson-D-meson interac- 
tions for all possible channels. In Section |2] the results in the 
physical world (i.e. at the physical values of the quark masses) 
are calculated to leading order (LO) and next-to-leading order 
(NLO) in the chiral expansion, as well as those obtained us- 
ing unitarized scattering amplitudes. In Section|3] we extrapo- 
late the results to unphysical quark masses using the unitarized 
scattering amplitudes. Some discussions and a brief summary 
is given in Section |4] Some useful isospin relations are rele- 
gated to AppendixlAl 
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2 Scattering lengths at the physical quark 
masses 

2.1 The chiral effective Lagrangians 

The leading order Lagrangian is just given by the kinetic and 
the mass term of the heavy mesons (chirally coupled to pi- 
ons) 012I13I14I 



QCD II16I17I . The he term is also dropped since it is suppressed 
by one order due to the commutator structure. We are therefore 
left with only two free, active parameters, namely and h^. It 
can be shown that the contributions of the and /13 terms to 
S-wave amplitudes differ only at 0{p/m£i). Therefore a vari- 
ation of /i3 within its natural bounds for a given provides an 
estimate for higher order contributions ITSl . 



(1) 



with D = (D", _D+, Df) collecting the charmed mesons with 

o 

the chiral limit mass Md- The co variant derivative is 



where 



U = exp 



F 



= U 



(2) 



(3) 



with F the Goldstone boson decay constant in the chiral limit 
(which we will identify with the pion decay constant in the 
following). The Goldstone bosons are collected in the matrix- 
valued field 



K- 



77+ K + 



(4) 



Since we are only interested in the scattering lengths with 
definite isospin, we only consider the strong interaction part 
of the NLO chiral Lagrangian describing the interactions of 
the pseudoscalar charmed mesons with the Goldstone bosons, 
which reads ITSl 



,(2) 



(5) 



where 



X+=X+ - 3 (x+> , 

The quark mass matrix is diagonal 

X = 2B ■ diag (m„, m^, 



(6) 



(7) 



in terms of i? — | (0|qq|0) The unknown coefficients = 
0, . . . , 6) are the pertinent LECs. At 0{p^), hi can be deter- 
mined from the mass differences among the D-mesons as hi — 
0.42 ifTJl . In the following, we choose to drop the terms with 
one more flavor trace in the NLO Lagrangian, i.e. the /iq, /i2 
and /14 terms, which are suppressed in the large A^c-limit of 



2.2 Numerical results for the scattering lengths 

We first consider the perturbative chiral expansion of the scat- 
tering amplitude, T = T^^^ ^T^"^^ + . . ., where the superscript 
denotes the chiral dimension. We work here to NLO, that is 
chiral dimension two. In that case, the scattering amplitudes 
are real and can be written as 



T{s,t,u) = T'^^\s,t,u)+T'^'^\s,t,u) 



4F2 



2Ci 
3^2 



hi 



2C: 



35 



F2 



i/35(s,i,M),(8) 



with 



h'iP2-Pi + hz {pi-P2P3.-Pi+Pi-PiP2-p?,) ■ (9) 



The coefficients in all the amplitudes for the Goldstone boson- 
_D-meson scattering are given in Table [T] The NLO S-wave 
elastic scattering amplitudes at threshold of the participating 
particles can be cast into the form 



T 1- 

^thi — p2 



2C 

C0M1M2 + -^hi + 2C35(/l3M| 



+2hzMlMl) 



(10) 



with Ml and M2 denoting the masses of the scattered heavy 
and light mesons, respectively. Note that from Table [T|one can 
see that the most attractive interaction occurs in the (S", /) = 
(1,0) DK channel, where S (I) denotes the total strangeness 
(isospin) of the two-meson system. This is the reason why 
in this channel the Z)*g(2317) was generated dynamically in 
many previous works II6I7I8I9I . 

The S-wave scattering length parameterizes the scattering 
amplitude at threshold 



ao 



87r(A/i + M2 



■Tthi ■ 



(11) 



Up to the order we are working, F can be replaced by the phys- 
ical pion decay constant = 92.4 MeV. We take the physi- 
cal masses for all the mesons, i.e., A/^ — 138 MeV, Mj^ = 
496 MeV, Af^ = 548 MeV, Md = 1867 MeV, and Md^, = 
1968 MeV 1 1 81 . The results for the S-wave scattering lengths 
at LO are given in the third column of Table |2] 

To work out the results to NLO, we take the dimensionless 
low-energy constant h'^ = h^M^a to be in the natural range 
of [—1, 1] as in Ref. ifTSl Pl CoiTespondinglv. h^ is determined 



' There is a typo in Ref. |151. There /15 was written as h^/M'^o 
while the correct one should be hsM^io ■ 
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Table 1. The coefficients in the scattering amphtudes. Here, S (I) denotes the total strangeness (isospin) of the two-meson system. 
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Table 2. The S-wave scattering lengths from calculations at LO and NLO (units are fm). The results using unitarized amplitudes are also given 
in the two columns denoted by UChPT and CUChPT, representing one-channel and coupled-channel unitarized chiral perturbation theory, 
respectively. 



from fitting to the mass of the Z3*q(2317) in the full calcula- 
tion, to be described below. This leads to /13 = —1.479 for 
/15 = 1, and /13 = 2.315 for /15 = —1. The results for the 
S-wave scattering lengths using these input parameters also for 
the perturbative calculation to NLO are given in the fourth col- 
umn of Table |2] with the uncertainties from the lack of knowl- 
edge of /13 and h'r,. We use these uncertainties as an estimate for 
possible higher order corrections. A comparison of the NLO re- 
sults with the LO ones shows a good convergence of the chiral 
expansion, especially for channels where only Z3-mesons and 
pions are involved. Even in the channels with kaons or etas, no 
dramatic NLO corrections are found. 

The above calculations are performed with perturbation the- 
ory up to a given order. However, a perturbative expansion to 
a finite order can never be reliable, if there is a bound state, 
a virtual state, or a resonance near by. All of these are non- 



perturbative phenomena, and the presence of such kind of state 
would modify the results from perturbation theory significantly. 
Therefore, we unitarize the amplitudes obtained from ChPT ac- 
cording to the method of Ref. |191 (called UChPT in the fol- 
lowing. An early review on this method is |20|). The unitarized 
amplitude is then given by the following resummation 

T{s)^V{s)[l-G{s)-V{s)]-\ (12) 

with V{s) the S-wave projection of the scattering amplitude 
T = r(i) + T(2) given in Eq. ©. Further, G(s) is the scalar 
two-meson loop integral lll9i21J regularized by a subtraction 
constant a{fi) with denoting the scale of the dimensional reg- 
ularization. The expression of the loop integral at threshold is 
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Table 3. Real parts, imaginary parts and Riemann sheet (RS) of the pole positions for the one-channel calculations for two parameter sets. All 
masses/energies are given in MeV. 
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Table 4. Positions of poles with the largest impact on physical observables from the coupled channel calculations for two parameter sets. 
Vii = denotes the results for a calculation where the diagonal interactions were switched off. Here the second (third) Riemann sheet for the 
(0, 1/2) case is defined by \m{qD-n)<Q, Im(gD77)>0, and Im(g£)^ j^')>0 (lm(go^)<0, Im(gD^)<0, and Im(gj3^^)>0), and that for the (1, 1) 
case is defined by Im(gD^-^)<0, and Im(gDx)>0 (Im(gDg^)<0, and Im(g_DK)<0). All masses/energies are given in MeV. 
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We will use the subtraction constant determined in Ref. [|T5|. 
i.e. a(l GeV) — —1.846. Then the scattering lengths with the 
unitarized amplitude for each channel can be obtained easily, 
and the results are given in the column denoted by UChPT in 
Table|2] For the channels with repulsive interactions, i.e. those 
with negative scattering lengths at LO and NLO in the conven- 
tion used here, the coiTection from the resummation is not dra- 
matic. It is several percent for the channels without kaons, and 
50% for the channels with kaons at most. For the channels with 
attractive interactions, i.e. those with positive scattering lengths 
at LO and NLO, the smallest correction is already 50% for the 
(0,1/2) Dt: channel. Both the scattering lengths of the (-1,0) 
DK and (0, 1/2) DgK channels change from about 0.3 fm to 
about 1 fm, and that of the (1, 0) DK even changes the sign. 
These dramatic changes can be understood on the basis of the 
singularity structure of the unitarized amplitudes. 

At this point some general remarks are in order Imagine 
some energy-independent potential that provides an attractive 
force between two particles. As the strength of the force is in- 
creased, eventually a virtual state appears in the unitarized scat- 
tering amplitude. A virtual state is a singularity below threshold 
on the real axis on the second sheet of the complex s-plane of 
the S-matrix. As the interaction strength is further increased 
the pole moves closer and closer to the threshold. At the same 



time the scattering length grows to become positive infinity 
when the pole hits the threshold. This situation is almost re- 
alized in nn scattering near threshold which shows an excep- 
tionally large scattering length of about 18 fm. If the strength 
of the interaction is increased even further, the pole jumps on 
the physical sheet — the state becomes a bound state. Corre- 
spondingly the scattering length changes sign and starts to in- 
crease from negative infinity as the pole moves away from the 
threshold — this is the case that applies to the deuteron chan- 
nel of nucleon-nucleon scattering, where the scattering length 
is about — 5 fm (Note that the sign conventions in nucleon- 
nucleon scattering are often different to those used here). These 
scenarios are discussed in various textbooks, see e.g. Refs. [22 23]. 
For energy-dependent interactions one may even start from a 
resonance, which corresponds to two poles in the second Rie- 
mann sheet, both located as mirror images with respect to re- 
flections on the real axis. This is illustrated as the points A and 
A' in Fig. [T] Then an increase of the interaction strength will 
let the magnitude of both the real part and the imaginary part 
decrease and will eventually lead to poles located below the 
elastic threshold, but still with non-vanishing imaginary parts 
(B and B' in the figure). A further increase in the strength of the 
potential might convert the two resonance poles into two virtual 
states (C and C in the figure). As the strength parameter gets 
increased even further, one pole moves towards the threshold 
to eventually become a bound state while the other one moves 
further away from the physical regime — this behaviour was 
e.g. reported in Ref. f24l for the case of the light scalar meson 
/o(600), where the interaction strength was varied by a varia- 
tion of the pion mass and was also observed in phenomenolog- 
ical studies 125 |. We checked that indeed the mentioned types 
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Fig. 1. Sketch of typical trajectories of S-matrix poles in the second Riemann sheet of the complex s-plane for energy-dependent potentials 
when some strength parameter is changed. See text for meaning of the labels. The dashed vertical line indicates the position of the elastic 
scattering threshold and the thick horizontal line the resulting unitarity cut. 



of poles appear by the unitarization. As can be read off Table[3j 
showing the pole positions for the single channel calculations, 
for /ig = +1 in the D^K channel there appears a pair of reso- 
nance poles below the DgK threshold. For the other parameter 
set displayed in the table (/ig = —1) those changed into two 
virtual states. For the full calculation with coupled channels for 
both parameter sets there are only resonance poles. In addition 
there are pairs of resonance poles in the second Riemann sheet 
for the channels (0, 1 /2) Dir and ( — 1,0) DK. The former one 
was found already previously in Ref. |7|. Note also that our 
findings for the scattering length in that channel are consistent 
with those of Ref. [1 IJ within the theoretical uncertainty. In the 
(1, 0) DK even a bound state appears. As a consequence the 
scattering length changes its sign. In Refs. 1261 in was shown 
that for bound states near a threshold one may write for the 
scattering length 




Applied to the i:)*o(2317) wehavee = Afi5+A/K-M£,.^(23i7) 
for the binding energy, /i — MdMk / {Mo + Mk) for the re- 
duced mass, and 1/(3, the range of forces, may be estimated via 
(3 ^ nip. The quantity Z, which can be identified with the wave 
function renormalization constant, is a measure of the molec- 
ular component of the state, with Z ^ \ ( Z — Q) for a pure 
elementary (molecular) state. Taking Z = 0, the above equa- 
tion gives a^ox^Dx = —1-05 fm, which is close to the value 
listed in Table |2l Thus, would this value be extracted from lat- 
tice simulations in the future, it would be a direct proof for the 
molecular nature of the £)*o(2317). 

For {S,I) = (0,1/2), (1,0) and (1,1), there are more 
channels with the same quantum numbers, and the unitarization 
should include also the coupled-channel effect (called CUChPT 
for simplicity). The off-diagonal interactions can induce an 
imaginary part to the scattering length for those channels cou- 
pling to a channel with a lower threshold. The results with 
coupled-channel unitarization are given in the column denoted 
by CUChPT in Table |2] From the amplitude coefficients for 
{S,I) = (0, 1/2) in Table [1] one finds the off-diagonal in- 
teractions for the DgK to Dir and Drj are large. As can be 
seen from the line labled by Vu = in Table |4] they alone 
can already produce poles. For the full calculation their ef- 
fect is to move the poles on the second Riemann sheet of the 



one-channel DgK — + DsK amplitude to a pole in the relative 
first Riemann sheet of the DgK channel (that is the Riemann 
sheet with positive imaginary part of the center-of-mass sys- 
tem (cms) three-momentum in this channel and negative cms 
three-momentum in the other two channels — for the full cal- 
culation this sheet is labled as III, for by convention the first 
sheet is related to the channel with the lightest threshold). The 
appearance of the pole on that physical sheet again leads to a 
change in sign of the real parts of the scattering lengths for 
both the Dt] and DsK channels. For {S, I) = (1, 1), the LO 
diagonal interactions vanish, and the NLO ones are repulsive. 
The off-diagonal interaction however produces a pole which 
is shown in the last line in Table |4] Thus the pole is purely 
a coupled-channel effect. Because of the pole the imaginary 
part of the scattering length in the (1,1) DK channel becomes 
much larger than its real part, a direct reflection of the coupled- 
channel effect. 

There are no experimental data for the scattering of Gold- 
stone bosons off Z3-mesons available. However, for four of the 
channels discussed above results of lattice gauge theory calcu- 
lations were recently published [|5( .These lattice data have to be 
taken with a grain of salt as they have been obtained for a single 
lattice spacing and volume, thus leaving room for sizeable sys- 
tematic discretization effects. The resulting scattering lengths 
are also shown in Table |2l To illustrate the level of agreement 
for the LO, NLO, and UChPT calculations with the lattice data 
the results for three of the channels are also shown in Fig. [2] 
(so far lattice data exist only for channels without channel cou- 
pling). For the fourth channel investigated on the lattice the 
scattering length is very small as it is for all the different calcu- 
lations we performed. Instead of this channel, in the figure we 
show the channel (—1,0) DK DK, which could be eas- 
ily investigated on the lattice with the same tools already em- 
ployed in Ref. 15 1. The first and important observation is that 
the NLO corrections are generally small — the series appears 
to be well behaved. In addition, the figure shows that for two 
of the channels, (0, f ) Dn Dtt and (2, i) DsK DsK, 
the agreement of the calculations with the lattice results is of 
similar quality for all three schemes, however, for the channel 
(—1,1) DK — > DK an agreement with the lattice results is 
achieved only after unitarization. This is a priori unexpected, 
since the interaction in this channel is repulsive. However, a 
repulsive interaction iterated to even orders gives an attractive 
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Fig. 2. Comparison of the scattering lengths in four channels with the lattice data. We give the results for LO, NLO and unitarized ChPT. 
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interaction and consequently by the unitarization the repulsion 
gets weakened — the scattering length gets smaller in mag- 
nitude. We checked that in our calculation the leading loop, 
formally a next-to-next-to-leading order contribution, already 
gives the bulk of the effect. Thus here the lattice results clearly 
show the trace of a higher order hadron-hadron interaction. As 
the first panel of the figure shows, we predict an even more dra- 
matic effect in the channel (—1,0) DK DK. As discussed 
above, here the interaction is attractive and the large change 
in the scattering length is the consequence of the appearance 
of a resonance state in the scattering matrix — clearly a non- 
perturbative phenomenon. It would therefore be of high theo- 
retical interest to have lattice results for this channel as well. 

Ideally lattice simulations would be available for all chan- 
nels calculated. As can be read off Table |2] in some channels 
the effect from the unitarization is very large while in others 
it is quite moderate. All this is calculated without any free pa- 
rameter. Therefore a comparison of the lattice results to those 
from our calculation would provide a very non-trivial test of 
the presence of non-perturbative hadron-hadron interactions. 

3 Chiral extrapolation 

The lattice results used here are calculated with finite lattice 
spacing and at unphysical masses of the u and d quarks — 
the masses of the strange and charm quarks are chosen at their 
physical values. In the lattice calculation of the charmed meson- 
light hadron interactions the lattice spacing is 6 = 0.12 fm. 
The mass of the s quark is 80 MeV, which is consistent with 
its physical mass, and four values are chosen for the masses 



of the u and d quarks which are 11 MeV, 16 MeV, 32 MeV 
and 48 MeV, respectively. Evidently, these are larger than the 
physical u {d) quark mass. To compare the lattice results with 
the physical quantities, one needs to do extrapolations. Here 
we only discuss the chiral extrapolation of the results from 
the unphysical quark masses to the physical quark mass. In 
Ref. |5|, the authors parameterize the scattering lengths using 
the first non-vanishing term in the chiral expansion. They used 
ci+C2Mt^ / to extrapolate the 7 = 3/2 Dtt scattering length, 
and ci -I- C2MI/FI to extrapolate the DsTr{K) and the / = 1 
DK^ scattering lengths. The values they got are shown in the 
last column of Table|2]and Fig. |2] 

The lattice results for the scattering lengths at unphysical 
quark masses are given at several chosen values of Mtt/Ftt, 
which were taken from Ref. f2E\. Therefore, we take the cor- 
responding values of unphysical M.^ from the same paper for 
the lattice data points shown in Fig. |3] Although the masses 
of the charmed mesons have been expanded to one loop or- 
der II29I . for our purpose we only need to work to 0{M^). 
From the Lagrangian given in Eq. (|5]l, modulo an overall un- 
measurable term proportional to (x+) which can be absorbed 
into the large-A^c suppressed ho term, one gets for the NLO 
correction to the mass squares of the charmed mesons 

SMl = AhiBm , SMl^ = ihiBrUs, (15) 

where rh = {rriu + md) /2. Since the lattice calculations |5 | use 
a value for iris which is close to the physical s quark mass, we 
will use the physical mass of the D^. From the first identity in 

^ In Ref. |5| this channel was called DK. However, the channel 
calculated is [cit][su], i.e. D°K+ (27]. 
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Fig. 3. Chiral extrapolation for the leading order results (dashed lines) and the full UChPT calculation (bands) compared with the lattice data. 



Eq. (fTSl l. the mass of the D meson up to NLO can be written as 

Ml 



Mr 



Ml 



hi- 



M, 



(16) 



using M^ = 2Brh. Using M|- = B{m + nis), one gets a 
similar expression for the kaon mass up to 0{M^) as 

M2 



M, 



K 



M 



K 



AMk 



(17) 



where M k denotes the mass of the kaon in the SU(2) chiral 

o 2 

limit rUu = mj, = 0, i.e. M ^ = Buris. With the above ex- 
pansion of Mo and Mk, one is ready to extrapolate the results 
given in Section |2] in the physical world to the world with un- 
physical quark masses. Assuming the subtraction constant to 
be A/jr-independent, which is true to NLO, the chiral extrap- 
olations of the results using UChPT for the same channels as 
in Fig.|2]are shown in Fig. [3] together with the lattice data. The 
description of the lattice data is rather good. The width of the 
band reflects the remaining freedom in the choice of the pa- 
rameters which is taken as an estimate of higher order effects. 
For comparison with the dashed line we also show the result to 
leading order. Although influencing the magnitude of the scat- 
tering lengths in all channels, only in the (0,3/2) Dtt chan- 
nel the higher order effects also changed the pion mass depen- 
dence, since only in this channel the changes in the light quark 
masses are not dwarfed by the large strange quark mass. 



As discussed before, in the (0, 1/2) D-k and the (1, 0) DK 
channels, poles get generated through the unitarization which 
govern the low-energy physics. Their presence is made appar- 
ent in Fig. m which gives the quark mass dependence of the 
scattering lengths for these channels predicted in CUChPT. To 
do the calculations the mass of the r/, which enters the coupled- 
channel calculations, was expanded as 



M^ = M^ 



Ml 

6M„ 



(18) 



The strength of the attractive interaction in the (0, 1/2) Dtt 
channel is measured by Mt^. When M^r is increased to some 
value, the interaction becomes strong enough to form a bound 
state. The structure shown in the left panel of Fig. |4]is a nice 
illustration of the physics discussed in Subsection 12.21 In the 
(1,0) DK channel, on the other hand, the strength of the inter- 
action is controlled by Mk, which is changed only moderately 
when the u, d quark masses get increased. As a result also the 
scattering length shows a weak A/^r dependence. 



4 Summary 

In this paper we calculated the scattering lengths for all the 
Goldstone boson-D-meson scattering channels using both chi- 
ral perturbation theory as well as a unitarized version of it. We 
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Fig. 4. Chiral extrapolation for the leading order results (dashed lines) and the full CUChPT calculation (bands) for the (0, 1/2) Dtt and the 
(1,0) DK channels. 



found that for some channels the unitarization produced new 
singularities in the S-matrix, in some cases as virtual states or 
resonances, in some cases as bound states. One of these dy- 
namically generated singularities can be identified with the ex- 
perimentally well established D*o(2317), that was identified 
as a hadronic molecule in various works 161718191 . However, 
we found more non-perturbative effects and demonstrated that 
a determination of scattering lengths on the lattice provides an 
ideal tool to investigate them. 
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A Isospin relations 

In the present work, we do not consider isospin violation. In 
this Appendix, we give the relations between the scattering am- 
plitudes in isospin basis and in particle basis which are useful 
in derivations. Before that, let us give the definition of Mandel- 
stam variables and the phase convention for the isospin eigen- 
states 

S={Pl+P2)^, t^{Pl-P3)^, W = (Pl -P4)^ , (A.l) 

and 

1 1' 



|^+) = -|l,+l), \K°) = - 
\D^) = - 



2' 2 



1 1 

2'^2 



(A.2) 



The two numbers in | - • •) on the rh.s. are /, I3. All the other 
states are defined with a positive sign. Then employing isospin 
symmetry and crossing symmetry, one can get the following 
isospin relations: 



V, 



DK^DK^ 
(-1,0) 



V, 



1^'"'^/'^ (S t U 



^(04/2) / ^ 
(04/2) 



V, 



(s, t, u 
(s, t, U 



(0,1/2) 



D,K^D 



is t u 



(1.0) 



DK^DK 



{s, t, u 



V, 



(1,0) 



V, 



1,0) 



V, 



D^-q^DK 
(1,1) 



{s, t, u 



V, 



(1,1) 



DK^DK 



(s, u 

y^''^/'' is t u 



= Vdok-~>dok- (s, t, u), 

-'2VD + K+^D+K+{'U,t,s) 

yD"K-^DOK- {Si ti u), 
= Vd + ti+~>D+tt+ (s, t, u), 

3 

= ■^VD+Tr + ->D+7T+ {u, t, s) 



D + ri~^D+ri 



{s,t,u), 



(A.3) 

(A.4) 
(A.5) 

(A.6) 
(A.7) 
(A.8) 
(A.9) 



= V3yo+^_^£,o^o(s,t, u), (A. 10) 

= ^D + K--»D«r,(^'*'")' (A. 11) 

= 2Vdo^-^£,o^- (m, t, s) 

{s,t,u), (A.12) 

= ~V2V^+.^^^0K+ {s, t, u)XA.U) 

= ^I? + 7rO^Ll + 7rO (A. 15) 

^VD + K+~>D + K+is,t,u), (A.16) 

^V2V^+j^-^jyo^oiu,t,s), (A. 17) 

= ^D + K+^D + K+('5>^''")' (A.18) 



where the superscripts mean {S, I) with S{I) representing the 
total strangeness (isospin) of the two-meson system. 
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